Constraining Sterile Neutrino Warm Dark Matter with Chandra Observations of the 
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We use the Chandra unresolved X-ray emission spectrum from a 12' — 28' (2.8-6.4 kpc) annular 
region of the Andromeda galaxy to constrain the radiative decay of sterile neutrino warm dark 
matter. By excising the most baryon-dominated, central 2.8 kpc of the galaxy, we reduce the 
uncertainties in our estimate of the dark matter mass within the field of view and improve the 
signal-to-noise ratio of prospective sterile neutrino decay signatures relative to hot gas and unresolved 
stellar emission. Our findings impose the most stringent limit on the sterile neutrino mass to date 
in the context of the Dodelson-Widrow model, rris < 2.2 keV (95% C.L.). Our results also constrain 
alternative sterile neutrino production scenarios at very small active-sterile neutrino mixing angles. 

PACS numbers: 95.35.-|-d, 13.35.Hb, 14.60. St, 14.60. Pq 
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I. INTRODUCTION 



In nearly all extensions of the Standard Model, the 
generation of neutrino masses leads to the introduction 
of sterile neutrinos, e.g., [Jill- Sterile neutrinos may also 
play the role of the dark matter, e.g., 0-13 and affect 
a host of interesting cosmological and ast rophys ical pro- 
cesses, including the production of baryon [Tsl - flTi [20j and 
lepton [3 asymmetries. Big Bang Nucleosynthesis ^ 
[2q | , the evol ution of the matter power spectrun i I2ll . [22 
reionization 
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neutrino oscillations [32|, |33| . pul- 
sar kicks and supernovae l4Qri2.]. The degree to 
which sterile neutrinos participate in these processes is 
sensitive to the strength of their interactions with Stan- 
dard Model particles and their abundance throughout 
cosmic history. The creation and exploration of models 
and production scenarios that characterize these proper- 
ties has been [3-7] and continues to be [8 12, 43-56] a 
very active area of research. See Ref. for a review of 
sterile neutrino properties. 

In addition to their rich phenomenology, sterile neu- 
trinos are also a readily testable dark matter candidate. 
To date, their properties have primarily been constrained 
in two ways: through X-ray searches for their radiative 
decays and via observed cosmological small-scale struc- 
ture, but the possibility of detection through atomic ion- 
ization and nuclear spin flip signatures has also recently 
been discussed [H^. Although the two primary meth- 
ods are generally considered separately, combined studies 
of both radiative and cosmological constraints have also 
been performed [i,!!^,!!^. 

The radiative decays of predominantly sterile neutrino 
mass eigenstates to predominantly active neutrino mass 
eigenstates and X-rays of energy = ms/2 can be de- 
tected by existing X-ray satellites [6|- Current radiative 
decay limits are based on observations of a lon g li st of 
sources, including the Cosmic X-ray background f6l','62"], 
nearby clusters such as Virgo [44] and Coma T3, 63], 
more distant clusters like A520 and A1835 [6J] and 
the bullet cluster [6^, nearby galaxies like Andromeda 



[66|,|67|, M33 [68|, Milky Way satellite galaxies like the 
Large Magellenic Cloud [69^, Ursa Minor [zO,[2l|, Draco 
[t^ ]. Wilman I [t!], and Seque I 74], as well as unre- 
solved emission from the Milky Way itself [zO, [75l - [78| . 
See Refs. fl3l,l66l.[76| and references therein for a more de- 
tailed summary and Refs. [7§,|8Q,] for more general studies 
of radiative limits. 

Although the original Dodelson-Widrow (DW) non- 
resonant sterile neutrino production scenario [3, @, 113] 
has nearly been excluded by radiative constraints alone, 
e.g, nis < 3.5 keV (95% C.L.) 66], and can account for at 
most 70% of the dark matter at 2 cr according to Ref. [5^ , 
recent Suzaku observations of Wilman I revealed a spec- 
tral feature at 2.5 keV that was consistent with the ra- 
diative decay of a 5 keV DW sterile neutrino 173 ] . T hese 
findings have been discussed further in Refs. [8l|,|82|, and 
we will return to them in Sect ions II VI and fVl of this paper. 

In addition to the DW scenario, many other models 
have been proposed that predict the proper relic dark 
matter density even at very small active-sterile neutrino 
mixing angles, e.g., [1, H, [IS IHj HI] • General phase space 
considerations have also been used to constrain sterile 
neutrino production scenarios, e.g., (sslls^l. 

As a complement to the upper bounds imposed by 
radiative decay constraints, observations of the cluster- 
ing of cosmological structures on the smallest scales can 
set lower bounds on the mass of the dark matter parti- 
cle. This is because sufficiently light dark matter parti- 
cles would have suppressed or even erased these smallest 
structures by easily propagating (free-streaming) out of 
the shallow gravitational potential wells in which they 
formed. 

Current cosmological lower limits on are based pri- 
marily on measurements of small scale clustering in the 
CMB, SDSS galaxy, and Lya forest flux power spec- 
tra and gravitational lensing by the smallest structures. 
However, the results of N-body simulations have also re- 
cently been used to set lower bounds on rus based on 
the requirement that the number of simulated satellite 
galaxies equals or exceeds the number of observed Milky 
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Way satellites [sBI]. 

Of all the cosmological constraints, the Lya limits are 
the most sensitive to the difficult-to-characterize, non- 
linear growth of baryonic and dark matter structures 
as well as the specifics of dark matter production, e.g., 
the momentum distribution of the dark matter particles. 
There have been significant discrepancies between the 
lower Hmits [13, [SGl - lssi lool . iQlj obtained using different 
data sets and hydrodynamical simulations, as discussed 
in, e.g., Lya limits are less restrictive for mod- 

els in which sterile neutrinos behave more like cold dark 
matter, e.g., p|. [TtI. [soj . 

Studies of gravitationally lensed galaxies and QSOs 
have provided comparable lower limits, e.g., > 5.2 
keV 01 and > 10 keV [o^l, respectively, based 
on the image distortions caused by the smallest struc- 
tures. These lensing constraints should improve signif- 
icantly as data from future submillilensing experiments 
become available [l^l- Lower bounds based on compar- 
isons of satellite galaxy counts to N-body simulations, 
e.g., rUs > 13.3 keV (in the DW scenario) ^85*], are also 
similar to the most restrictive Lya constraints. 

In this paper, we consider the radiative decay limits 
imposed by Chandra [9^ observations of the Andromeda 
galaxy (M31). We choose to focus on Andromeda be- 
cause of 1) its close proximity, 2) its substantial and well- 
studied dark matter distribution, and 3) its intrinsically 
low level of X-ray emission. These properties tend to 
maximize the prospective sterile neutrino decay signal (1 
and 2) and minimize noise, i.e., astrophysical background 
(3). 

The Chandra data set we use in this study also has 
several advantages over the XMM-Newton observations 
of Andromeda 96] that were used to constrain the prop- 
erties of sterile neutrinos in Ref. [66*1 (W06). First, com- 
pared to XMM-Newton EPIC, the Chandra ACIS detec- 
tor has a significantly lower and much more stable instru- 
mental background, making it particularly well-suited for 
extracting low surface brightness, extended X-ray emis- 
sion. Second, using Chandra's superior angular resolu- 
tion, we are able to remove the emission from a larger 
number of resolved. X-ray point sources than would be 
possible with XMM-Newton, thereby further reducing as- 
trophysical background noise, while excising very little 
dark matter from the small excluded regions. Third, the 
field of view (FOV) associated with the Chandra spec- 
trum we use in this paper is over 25 times the area of the 
XMM-Newton FOV studied in W06 (a 12 arcminute to 28 
arcminute annulus vs. a circle of 5 arcminute radius) and 
therefore probes a significantly larger dark matter mass 
and prospective decay signal. The dark matter mass 
within this FOV is not only significantly larger but also 
subject to less uncertainty both because we exclude the 
central region of the galaxy, where the dark matter den- 
sity profile is least well constrained, and because we make 
use of an updated model of the Andromeda dark matter 
distribution based on more recent kinematic data and 
theoretical considerations. The results of this study are 



also more accurate because we use the Chandra ACIS-I 
effective area, rather than an average flux-to-count ratio, 
to convert calculated sterile neutrino decay signals into 
spectral features. For Majorana sterile neutrinos, the re- 
sulting unresolved emission spectrum requires ms < 2.2 
keV (95% C.L.) to avoid more than doubling the ob- 
served signal in bins of energy E > 1.1 keV. Although 
the mass-mixing exclusion region we generate with this 
new data set is not vastly more restrictive than the ex- 
clusion region presented in W06 (see Fig. 4 below) , it is 
more robust for the reasons cited above, and as we argue 
throughout the paper, it will be difficult to improve upon 
the constraints we set here with the current generation 
of X-ray detectors. 

In Sec. ini we discuss the basics of the DW scenario 
as well as some more recent sterile neutrino models. In 
Sec. mil we discuss the Chandra observations of An- 
dromeda, how we generate the unresolved X-ray spec- 
trum we use for our constraints, and how we calculate the 
dark matter mass contained within the Chandra FOV. In 
Sec. llVi we present our analysis of the unresolved spec- 
trum, discuss the limits we are able to impose on the ster- 
ile neutrino mass and mixing, and compare our results to 
those of previous studies. In Sec. |V] we summarize our 
findings and discuss the outlook for future sterile neu- 
trino constraints. Throughout the paper, we assume a 
flat cosmology with fibaryon = 0.04, r^vvDM = rig = 0.24, 
Qa = 0.72, and h = i/o/100 km s"^ Mpc^i = 0.72. 



II. THE STERILE NEUTRINO WARM DARK 
MATTER MODEL 



The radiative decay rate for Majorana sterile neutrinos 
is MM 



Ts ~ 1.36 X IQ-^^s-i 



sin^ 261 



keV 



(1) 



The line flux at E^y^s = rns/2 resulting from the decay 
of the fraction of sterile neutrinos that are within the 
detector field of view (FOV), A^f°v ^ (M™//^^) a 
halo at a distance D is given by |3j 



$,,s(sin" 29) 



2 on^ _ ^7.siV, i . ^ X lO-i^erg cm-2s-i 



AfFOV 



lO"M0y VMpc 



D 



sin^ 29\ / m, \5 



10 



-10 



where sin^ 29 characterizes the active-sterile neutrino 
mixing. Dirac sterile neutrinos would produce only half 
the flux given by Eqn. © [97, 98]. To facilitate compar- 
isons between our work and other results, we will adopt 
Majorana sterile neutrinos, which have been more com- 
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monly assumed in recent studies, e.g., |67l l76l l7a| . 

For a QCD phase-transition temperature of Tqcd = 
170 MeV and a lepton asymmetry of L ~ nharyon/n^ ~ 
10""'^'^, the sterile neutrino density-production relation- 
ship for this model is ^] 



TUs — 55.5 keV 



sin226l 



-0.615 



0.24 



0.5 



(3) 



Combining Eqns. © and ([3]) yields an expression for the 
line flux that is independent of the mixing angle: 



7.0 X 10"^^erg cm^^s ^ 



VioiiMgy VMpc 



0.24 



0.813 



V keV/ 



3.374 



(4) 



We note that the density-production relationship pre- 
sented in Ref. [i^ agrees with Eqn. ^ for sterile neu- 
trino masses ranging from 1 keV ^ nig S 10 keV. The 
sterile neutrino mass limits we determine in Sec . IIVI are 
therefore valid for the relationships in both Refs. [4j] and 
[48l |. It would also be straightforward to re-evaluate our 
limits on assuming different density-production rela- 
tionships, such as the three resonant production models 
associated with large L [l^ or the Shi-Fuller model f3| 
as calculated in Ref. [11] , all of which are shown in Fig. |3] 
below. 



III. PROPERTIES OF ANDROMEDA 




FIG. 1: Here we show the raw counts associated with the 7 
Chandra ACIS-I exposure regions (some of which are overlap- 
ping) within a 12' — 28' annulus about the center of the An- 
dromeda galaxy. Variations in the brightness among these 7 
regions are mainly due to differences in exposure times (which 
vary from 5 ks to 20 ks) rather than intrinsic variation in the 
emission (see text). We have also included an ACIS-I image 
of the central 12' to illustrate the extent of the X-ray emission 
from hot gas and point sources that we have eliminated by 
excising the bulge of the galaxy from our FOV. 



A. X-ray Data 

To detect or place the most restrictive constraints on 
radiatively decaying dark matter, it is critical to mini- 
mize the X-ray background from baryonic sources, such 
as diffuse hot gas and X -ray binaries. In M 31, b oth the 
diffuse hot gas [H, Il00j | and X-ray binaries |l0ll | exhibit 
a relatively steep radial distribution toward the galactic 
center, suggesting that we can optimize the dark mat- 
ter signal to baryonic noise ratio by avoiding the inner 
regions of the galaxy. Although most existing Chandra 
observations of M3I have been aimed toward the inner 
bulge, we did find seven observations that were aimed 
toward the outer bulge/disk regions and were therefore 
more well-suited for our purposes. All seven observations 
used the Advanced CCD Imaging Spectrometer (ACIS) 
I-array as the primary detector, with a 17'xl7' field of 
view (FOV) . The positions of these observations relative 
to an annulus of inner radius 12 arcminutes and outer ra- 
dius 28 arcminutes ( 12' - 28' ~ 2.8 - 6.4 kpc at Dmsi ^ 
0.78 ± 0.02 Mpc 3M [lOl) are shown in Figure [H Five 
of the seven observations (ObsIDs 1576, 1584, 2899, 2901 



The decay rate equation given in Ref. [76l | is for Majorana rather 
than Dirac sterile neutrinos. 



and 2902) have effective exposure times of ~5 ks, while 
the other two observations have longer exposure times 
of 8 ks (ObsID 7067) and 20 ks (ObsID 11256), respec- 
tively. Each 17'xl7' ACIS-I region covers 11.4% of the 
12'— 28' annulus, and the non-overlapping parts of these 7 
pointings collectively cover ^ 35%. We address sampling 
bias later in this section, pointing out that the co-added 
spectrum from the selected observations provides, if any- 
thing, an overestimate of the X-ray emission from the 
entire 12' — 28' annulus. 

We obtained and re-processed the archival data using 
CIAO version 4.2 and the corresponding calibration files. 
We followed the standard procedure of data reduction, 
e.g., as described in Ref. [99|. Briefiy, for each observa- 
tion, we (i) filtered time intervals of high particle back- 
ground, (u) detected and removed discrete sources, (in) 
extracted a spectrum for the unresolved X-ray emission in 
the 0.5-8 keV range and generated corresponding instru- 
mental response (rmf and arf) files, and (iv) extracted 
a fiducial instrumental background spectrum, using the 
"stow background data" that has been calibrated with 
the 10-12 keV count rate. In a final step, we co-added the 
seven spectra and generated exposure- weighted response 
files, using the FTOOL addspec. Such a procedure is ap- 
propriate, since the spectra were extracted from nearly 
identical detector regions. 
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Since the effective area Aes{E) of ACTS is known to 
decrease gradually with time, it is important to note that 
these seven observations were taken over a time span of 
9 years. By using a weighted Acs{E) associated with the 
co-added spectrum in our analysis below, we are able to 
gauge the mean sensitivity of the detector over this time 
interval. Additionally, we note that the degradation of 
the detector mainly affects energies below 1 keV, and our 
data set has very little constraining power at these low 
energies, as we will show. 

It is also noteworthy that the observed flux varies mod- 
estly among the seven observations, within a factor of 1.5 
of the mean values associated with the co-added spec- 
trum. This is a natural result, since the observations 
sampled different regions of the outer bulge and the disk. 
We note, however, that the sampled regions mostly cov- 
ered the southeastern half of the disk - where the unre- 
solved X-ray emission is higher than that from the north- 
easter side, primarily due to the tilt of M31's galactic 
plane [99|]. This sampling bias means that the co-added 
spectrum we use represents an upper limit on the mean 
flux within the 12' — 28' annulus - precisely what we need 
to set a robust upper limit on the mass of the sterile 
neutrino. Finally, we note that a correction factor of 
(0.114)^-'^ has been applied to the data shown in Fig. [2] 
to account for the fact that the co-added spectrum from 
the 7 archived observations represents 11.4% of the pro- 
jected area of the 12' — 28' annulus. 



B. Dark Matter Enclosed within the Chandra Field 
of View 

The dark matter halo of M31 has been studied exten- 
sively. Klypin, Zhao and Somerville (KZS) jl04| . for in- 
stance, used a variety of kinematic data to determine its 
dark matter density profile, /9dm,M31- Using more recent 
data and a more accurate baryonic mass profile, Seigar, 
Barth and Bullock (SBB) [105] have updated the pro- 
file found by KZS. When we integrate the SBB model of 
Pdm,M31 over the volume, Vpov.Msi, associated with the 
projection of the Chandra FOV along the line of sight, 
we find that the 12' — 28' annulus contains: 



^FOV 
■^DM.MSl 



(0.8 ±0.1 



X lO"M0Mpc"^; 



where 



A^FOV 

-""DM.MSl 



^FOV 



(0.49 ±0.05) X lO^Mg 



pDM{\r-D\)dVi 



FOV 



(5) 



(6) 



is the dark matter column density within VpoV: and the 
dark matter mass within V^fov is defined by M^'^ ~ 
D^sgov^ (ggg for further details). 

In addition to being based on more contemporary 
data and theoretical considerations, the SBB value for 
^dmVsi is also more conservative: less than 85% of 
the KZS estimate for this annulus. We further note 



that for our FOV, CDM models predict more conserva- 
tive Sp^jyjg^ values than cored, WDM models. Because 
WDM models are less centrally concentrated than CDM 
models, they are constrained to predict higher densities 
at larger radii in order to reach agreement with the to- 
tal mass of a given halo. For instance, within the radial 
limits of our FOV (2.8-6.4 kpc), the dark matter density 
given by the Burkert profile for M31 [106] with the com- 
monly adopted parameter values (as in, e.g., Ref. [67j ) 
is 17% - 29% higher than that of the KZS profile and 
18% - 41% higher than that of the SBB profile. There- 
fore, despite the fact we are constraining a WDM can- 
didate in this paper, we chose the SBB profile over a 
WDM profile in the interest of providing a conservative 
dark matter mass estimate. 

Because of Chandra's excellent angular resolution, we 
lose less than 5% of S^^j^g^ due to the excision of 
point sources. However, to be conservative, we use only 
0.95 Sd^^m31 for ^DM in Eqns. and dH) to determine 
the prospective sterile neutrino signals in the figures be- 
low, and we ignore the contribution from the fraction of 
the Milky Way dark matter halo within the FOV. 



IV. CONSTRAINING STERILE NEUTRINO 
DECAYS 



In Fig. [21 we show the unresolved emission spectrum 
from the inner 12' — 28' of Andromeda. To calculate 
the Vs decay signals, we assume that sterile neutrinos 
comprise all of the dark matter, i.e., fig — SIdm = 0.24, 
and evaluate Eqn. ^ based on 95% of the Sj^^j^g^ value 
given in Eqn. ([S]). We then convert the sterile neutrino 
decay fluxes to the same units as those of the measured 
spectrum (Counts/sec/keV) as follows: 



E. 



6.7 X 10^2 Counts/sec/keV 



Acff (E^^s) 



VFOV 



10"MqMpc-2 J VO.24 



100 cm2 



a \°-^" / m, \ 1-374 



where 1 eig/E^^s — 1-6 x 10^(i?^_s/keV)~^ gives the num- 
ber of Counts associated with 1 erg at energy Ej^s and 
Ach{E~^^s) and AE are the effective area^ and spectral en- 
ergy resolution of the Chandra ACIS-I detector, respec- 
tively. To realistically simulate detected sterile neutrino 
decay "line" fluxes in Fig. [51 we use a Gaussian centered 
at E^^s = ms/2 with a FWHM of AE = E^^^fi^^, a con- 
servative estimate of the energy resolution of ACIS-I in 
imaging mode"^. 



^ http://cxc.harvard.edu/proposer/POG/html/ACIS.html 
^ http:/ /heasarc. nasa.gov/docs/cxo/cxo.html 
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TABLE I: The Table shows the distance to Andromeda, the 
angular range, ASpov, of the annular field of view (FOV) 
probed by the Chandra observations, the dark matter col- 
umn density, Ej^^, enclosed within the FOV (Eqn. [SJ, and 
the (95% C.L.) upper bounds on nis for Dirac and Majorana 
sterile neutrinos. 



Galaxy Name 



Andromeda (M31) 



Distance (Mpc) 



0.78 ±0.02 



ASpov (arcminutes) 



f 2' - 28' 



EE?,^/10"MeMpc-2 



0.8 ±0.08 



(keV) (95% C.L.) 



2.4 (Di 



mf (keV) (95% C.L.; 



2.2 (Majorana) 



A. DW Sterile Neutrino Mass Limits 

To determine the mass limit imposed by the unresolved 
emission spectrum of Andromeda, we find the first bin 
of energy E^^^ = ^sl"^ for which the sterile neutrino 
decay signal (Eqn. [71) at least doubles the amplitude of 
the measured spectrum in that bin, AJ^: 



7,s 



(a) > A^. 



(8) 



The resulting limits: to^ > 2.2 keV (Majorana) and 
> 2.4 keV (Dirac), which are shown in Fig. 2 and 
Table I, are much more significant than the 95% C.L. de- 
fined by the (1 a) Poisson error bars on the measured 
points. (For this data set, Eqn. ^ is equivalent to 
4 cr). However, because we do not have a precise un- 
derstanding of the features of the unresolved emission 
spectrum, which originate from some combination of hot 
gas, unresolved stellar sources, sky and detector back- 
grounds, and possibly sterile neutrino decay, we choose a 
limiting criterion that would remain robust even to 100% 
level fiuctuations in the data. (See W06 [i^ for a more 
detailed discussion). 



B. Examining Possible DW Sterile Neutrino Decay 
Signatures 

There are a large number of atomic emission lines 
at energies just below our exclusion limit, i.e., E^^s'^ 1 
keV. Differentiating between these features and anoma- 
lous lines will be difficult if not impossible given the spec- 
tral resolution of current detectors. The diminishing am- 
plitude of sterile neutrino decay signatures at lower en- 



oc (Eqn. [T]) further exacerbates this 



problem. Even with a data set such as ours, with a fa- 
vorably high ratio of FOV dark matter to astrophysical 
background, prospective sterile neutrino signals become 
comparable to or dwarfed by atomic line features at en- 
ergies just below 1.1 keV. 

To illustrate this point, we examined whether or not 
any of the atomic emission lines below 1.1 keV in Fig. 2 
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FIG. 2: Here we show the Chandra unresolved X-ray spectrum 
emitted from a 12' — 28' annular region about the center of the 
Andromeda galaxy (solid) and the first statistically significant 
Vs decay peaks (dashed) at E-y^s = ms,iim/2 = 1.1 keV (blue) 
and 1.2 keV (red), which exclude > 2.2 keV (Majorana) 
and m° > 2.4 keV (Dirac), respectively (95% C.L.; Eqn.[8|. 
As a gauge of the statistical significance of our limits, we 
have included the (1 a) Poisson error bars on the unresolved 
emission spectrum data points. As a point of comparison 
to the potential detection discussed in Ref. 73], we also show 
the decay signature that would be produced if the dark matter 
within the 12' ~ 28' region of Andromeda were composed of 
5 keV sterile neutrinos (green, dot-dashed), a possibility that 
is strongly excluded by the data. 



were statistically consistent with sterile neutrino decay. 
We found that this was the case for the Ne IX peak at 
1.07 keV. In particular, when we assumed that the Ne 
IX peak was a sterile neutrino decay line superposed on 
a continuum background characterized by the mean fiux 
values of the bins to the left and right of the peak energy, 
we found that it was best fit by the decay signature of a 
2.13 keV Majorana sterile neutrino (see Fig. [3]). 

This exercise underscores both the spectral resolution 
and amplitude issues. If Andromeda's dark matter halo 
were actually composed of 2.13 keV sterile neutrinos, 
they would produce a decay signature at 1.065 keV, but 
this is clearly indistinguishable from the 1.07 keV Ne IX 
feature with the detectors aboard Chandra. If the sterile 
neutrinos were of even lower mass, the amplitude of their 
decay signatures would be dwarfed among the thicket 
of sub-keV atomic lines, and higher spectral resolution 
would become indispensible if we were to have any hope 
of finding them. One of the few options available in the 
absence of such a dvan ces is to consider anomalous line 
ratios, as in Ref. (l07j |. but doing so requires extremely 
precise knowledge of chemical abundances, plasma tern- 



6 



> 

CD 
J4 



o 

CD 



-i-> 

:=! 
o 



1.02 




100 



E [keV] 



FIG. 3: Here, to illustrate how difficult it is to distinguish 
between atomic and anomalous line features with current de- 
tectors, we show the statistical consistency of the 1.07 keV 
Ne IX emission peak {Chandra data points in black) and the 
decay signature of a 2.13 keV Majorana sterile neutrino (blue, 
shaded region). The light blue band shows the 1 a uncertainty 
range associated with the Chandra data. 



perature, etc. of the region(s) of the target galaxy being 
examined. To make progress without such complications, 
particularly at E^^s"^ 1 keV, a new generation of much 
higher spectral resolution detectors is required, as we dis- 
cuss further in the Sec. IVl 



C. Exclusion Regions in the Mass-Mixing Plane 

To determine the region of the Ws — s.it?20 (mass- 
mixing) plane (Fig.Hl) that is excluded by the unresolved 
X-ray spectrum of Andromeda, we convert Eqn. ^ to 
Counts/sec/keV: 



dN, 



7,s 



dEj^sdt 



sin^26l = 



$,,s(sin2 20)\ fA,s{E, 



E. 



7:S 



9.8 X 10"^ Counts/sec/keV 



Acff (E-y^s) 
100 cm2 



VFOV 



101iMqMpc-2 



and adopt the analog of Eqn. 
rion: 



sin^ 20 
as our exclusion crite- 



dN, 



7,s 



dEjgdt 



3in^ 261) > A J". 



(10) 



Just as we found two mass limits, we also derived two 
exclusion regions for Dirac and Majorana sterile neutri- 
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FIG. 4: Here we present constraints on rus as a function of 
mixing angle, sin^2S, assuming that all dark matter is com- 
prised of sterile neutrinos (ils = 0.24). For L ~ 10~^°, the 
thick, solid line corresponds to Q,^^ — 0.24 ± 0.04 in the 
Dodelson-Widrow (DW) scenario (Eqn. [3]), while the region 
to the right corresponds to > 0.28. Three density- 

production relationships associated with ilg = 0.3 and (left 
to right) L = 0.1, 0.01, and 0.003 are also shown (dotted) 
, as is the Shi- Fuller density-production relationship com- 
puted in Ref. [5^ (dashed). The three previous radiative de- 
cay upper limits (all 95% C.L.) are based on Integral mea- 
surements of the unresolved X-ray emission from the Milky 
Way halo '77, 78], HEAO-1 and XMM-Newton observations of 
the Cosmic X-ray Background (CXB) [6l|| , and the most strin- 
gent constraints from the many limits imposed by nearby 
galaxies and clusters [H, H, HI, HI, IH, [t^, [zgl . The magenta 
line shows the recalculated boundary of this exclusion region 
for Majorana sterile neutrinos (see text), allowing for direct 
comparison between our results and those of Ref. [6^. The 
most restrictive radiative decay limits, from the present work 
(also 95% C.L.), are based on Chandra observations of the 
Andromeda galaxy. 



nos. The most restrictive region, which was determined 
by comparing the Chandra unresolved X-ray spectrum 
to the Majorana sterile neutrino decay flux, is shown in 
Fig. m The "indentation" of our exclusion region at the 
highest masses comes about because the effective area of 
the ACIS-I detector falls even more steeply than the spec- 
tral data at the highest photon energies. As discussed in 
the conclusion, a new instrument with a much larger ef- 
fective area and superior spectral energy resolution will 
be required to dramatically improve upon the radiative 
decay constraints presented here. 

In addition to our new Andromeda bounds, (the dis- 
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tinct parts of) three previously determined radiative de- 
cay exclusion regions are also shown in Fig. 21 The upper 
exclusion region is based on Integral measurements of the 
unresolved X-ray emission from the Milky Way halo [t^] ■ 
We note that these results were corroborated by a very 
similar later study [ill . The second exclusion region was 
derived by analyzing HEAO-1 and XMM-Newton obser- 
vations of the Cosmic X-ray Background (CXB) j3i|. The 
third region represents the most stringent constraints (66| 
(W06) from the many limits im p ose d by nearby galaxies 
and clusters 0, 113, 113, IH, [eg [t^, [zg. We note that 
if we were to recalculate the exclusion region derived in 
W06 [6^ for Majorana rather than Dirac sterile neutri- 
nos, the boundary of the exclusion region would remain 
the same (the magenta line shown in Fig. 4). This is 
because we would not only need to multiply $x,s by a 
factor of 2, we would also have to change the overly op- 
timistic estimate of the spectral energy resolution used 
in that paper {AE — E/30) to the more realistic value 
used in this analysis {AE = E/15). Since ^'^^ cx 

(Eqn. 9), the factors of 2 offset. 

Boyarsky et al. also conducted a very thorough 
analysis of an annular (XMM-A''euitoTi-observed) region of 
Andromeda (5' — 13'). Unfortunately, because of XMM- 
Newton's poorer spatial resolution, they were forced to 
excise almost one fourth of the field of view 23%) to 
remove the emission from point sources, thereby sacri- 
ficing potentially signal-producing dark matter to reduce 
the astrophysical background. As a result, the unresolved 
emission spectrum they generated probed a much smaller 
dark matter mass than the Chandra unresolved X-ray 
spectrum we consider here, and their limits are corre- 
spondingly less restrictive (rn^ < 4 keV). 

V. CONCLUSIONS 

In this paper, we used the Chandra unresolved X-ray 
spectrum of the Andromeda galaxy (M31) to improve 
the radiative decay constraints on sterile neutrino warm 
dark matter. Assuming either the model described in 
Refs. d, H, m, nil or Ref. [Ill, our analysis requires 
nis < 2.4 keV (95% C.L.), in the least restrictive case 
(Dirac sterile neutrinos). Because many recent papers, 
e.g., [13, [zl, [tII have assumed Majorana sterile neutri- 
nos, we quote the Majorana limit: 

m,s < 2.2 keV (Majorana; 95% C.L.) (11) 

as our lower bound for the DW scenario, so that compar- 
isons to other studies can be made on equal footing. 

Most of the currently available sterile neutrino mass- 
mixing parameter spac e remains open only for scenarios 

* |http://ixo.gsfc.nasa.gov| 



in which fls ^ 0.3 can be generated at very small mixing. 
Sterile neutrinos that are, e.g., resonantly produced in 
the presence of a large lepton asymmetry {L » 10"^"), 
e.g., 0, H, [H, [11] or created via Higgs decays [io^ are 
still viable. However, the combination of our X-ray con- 
straints with those from nearby sources are also able 
to partially restrict the L — 0.003 resonant production 
model [13] and the Shi-Fuller model [5^] for sterile neu- 
trino masses in the 10 to 24 keV range (Fig. 4), and the 
Integral measurements of the Milky Way halo [tJ [7§| ex- 
clude all the alternative production scenarios shown in 
Fig. m above = 40 keV. 

Even without cosmological small scale structure 
bounds, the DW scenario remains viable only between 
the Tremaine-Gunn bound [lOSj and our limit from this 
work, 0.4 keV < rus < 2.2 keV, which interestingly falls 
within the range of dark matter particle masses that 
best e xplain s the core of the Fornax Dwarf Spheroidal 
galaxy [l09| . This result underscores the need to continue 
to carefully and independently pursue all possible con- 
straints on sterile neutrino properties. However, as noted 
above, because of the decreasing sterile neutrino signal at 
lower rUs values and the large number of atomic emission 
features at energies ^ 1 keV, it will be difficult to improve 
upon the 2.2 keV limit we have presented here with ex- 
isting X-ray detectors. One of the most promising routes 
toward improved radiative constraints includes the use of 
much higher spectral resolution instruments than those 
currently available, as discussed, e.g., in Refs. fidi. IllOl - 
Ill2l |. The International X-ray Observatory (IXO^), for 
instance, which is scheduled for launch in 2021, will have 
a FOV comparable to Chandra detectors (~ 18' for pho- 
ton energies ranging from 0.1-15 keV), but ^ 100 times 
their effective area, ~ 10 times their spectral resolution, 
and ~ 10 times lower instrumental background. Based 
on these specifications, a megasecond observation of An- 
dromeda with IXO will have the capacity to fully test the 
4 alternative sterile neutrino production scenarios shown 
in Fig. |4] over the entire range of m ass-in ixing parameter 
space for which they remain viable |llll |. 
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